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I We discuss a system comprising two nonlinear (Kerr-like) oscillators coupled mutually by a 

nonlinear interaction. The system is excited by an external coherent field that is resonant to the 
frequency of one of the oscillators. We show that the coupler evolution can be closed within a 
. , finite set of n-photon states, analogously as in the nonlinear quantum scissors model. Moreover, 

I for this type of evolution our system can be treated as a Bell-like states generator. Thanks to the 

' nonlinear nature of both: oscillators and their internal coupling, these states can be generated even 

o ■ 

\^ . if the system exhibits its energy dissipating nature, contrary to systems with linear couplings. 
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I. INTRODUCTION 



In recent years the interest in the possibihty of generation and manipulation of non- 
classical states has been growing. It should be stressed that such a possibility is very 
important for the implementation of models in the quantum information theory. In fact, of 
special interest are the systems, whose dynamics can be closed within a finite set of n-photon 
states, due to the potential possibilities of implementations of the models discussed to the 
quantum computational systems. The systems leadingto the finite-dimensional state gen- 
eration are referred to as linear 0,|2| or non-linear ^,1^1 quantum scissors. Also, nonlinear 
couplers (NC), i.e. the systems involving two Kerr nonlinear oscillators with linear interac- 
tions between them, coupled with an external field have been considered in that context P| . 
The idea of the NC was introduced by Jensen ^ and Maier ^land has been developed and 
investigated in many aspects (for a comprehensive review see jsf). For instance, NC systems 
can exhibit such phenomena as: self-modulation, self-switching, self-trapping chaotic 
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systems synchronisation |l^ and others. Moreover, NC can ex. 
features such as sub-Poissonian or squeezed light generations 
interesting from our point of view, these models can also lead to generation of the finite 
dimensional and maximally entangled states. As we shall show in this paper, the Bell-like 
states can also be produced by the NC systems. However, one should realize that the sys- 
tems containing the Kerr nonlinearities are very sensitive to the damping processes |15l and 
hence, the losses of various nature are able to destroy generation of the states discussed 
almost completely. Nevertheless, we shall show that thanks to the nonlinear nature of the 
internal coupling involved, the influence of the damping is not so pronounced as for the 
systems with linear couplings jlG^ . 



II. THE MODEL AND SOLUTIONS 



The model of a nonlinear coupler considered here is built of two nonlinear oscillators 
characterized by Kerr nonlinearities Xa and xt, corresponding to the two cavity field-modes, 
labelled by a and b respectively. These oscillators are located inside a high-Q cavity and are 
coupled together by the interaction of the nonlinear character. Additionally, the field mode 
corresponding to the one of the Kerr oscillators (mode a) is linearly coupled to an external 
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classical field. Thus, the system discussed can be described by the effective Hamiltonian (in 
the interaction picture): 

H = HnL + Hint + Hext , (1) 

where 

Hnl = ^{a^ra' + fibYb' , (2a) 

H,^t = e{dYb' + ^*{b^fa' , (2b) 

Hext = cta'' + a*a . (2c) 

a(a^) and 6(6^), appearing above, are the usual bosonic annihilation (creation) operators, 
corresponding to the a and b modes of the oscillators, respectively. The parameter e describes 
the strength of the internal coupling between the two oscillators, whereas a is the strength 
of the linear coupling between the external classical field and the field of the mode a inside 
the cavity. 

Since this section is devoted to the ideal situation, i.e. the model without damping 
processes, we shall describe the system evolution in terms of the time-dependent wave- 
function (a more realistic situation of a damped coupler will be also considered and discussed 
in the fifth section). It can be expressed in the n-photon Fock basis as: 

oo 

I^W) = ^"''"(^^ l'^^'^ \^^b ' (3) 

n,m=0 

where Cn,m{t) are the complex probability amplitudes of finding the system in the n-photon 
state for mode a and the m-photon state for mode b. 

Since our model involves the external pumping by the classical coherent field, the total 
energy of the system considered is not conserved. Therefore, one could expect that with 
increasing number of photons the Fock states will be involved in the system dynamics. How- 
ever, the assumption of the weak external pumping and constant amplitudes of the couplings 
enables us to use the method of the non-linear quantum scissors extensively discussed in 
0, In consequence, we are able to truncate the wave function (jH)) to the wave function 
describing only the evolution of the resonant states. The form of Hnl implies that only 
|0)^ |2)^, |1)^ |2)^ and |2)^ |0)^ states should be considered. Thus, the wave function takes 
the form: 

\m)cut = C2,o(t) |2), |0), + c,,2{t) |1), |2), + coAt) |0)„ |2), , (4) 



where we have restricted our considerations to a finite set of the states. 

Applying the Schrodinger equation we are able to write the set of equations of motion 
for three probability amplitudes in the closed form (we use units of h = 1): 

■d 

i-r:C2,o = 2eco,2 , 
at 

d 
dt 

^^co,2 = 2e* C2,o + a* ci,2 • (5) 

Assuming that for the time t = we have two photons in mode a and no photons in mode 
b inside the cavity i.e. lipit = 0))^^^ = |2)^ |0)^, we can easily solve the set of equations (jS)) 
showing that the solutions have the form: 

a'^ + Ae^ cos(nt) 



Cl,2(t) 

co,2(t) = -2t \ ' , (6) 



1]2 

ea (cos(i7t) 
e sin(f2t) 



where we have introduced the effective frequency Vt = {a^ + 4e^)^/^. 

To validate our analytical results we can also perform numerical calculations to find 
the probability amplitudes using the method described in ^]. Thus, we construct the 
unitary evolution operator JJ = exp (^—iHt^ using Hamiltonian ^ and apply it to the initial 
system state \ip(t = 0)). For the situation discussed here, we obtain the wave function in 
the following form: 

\m) = m\2)jo), . (7) 

What is important for the numerical calculations performed here, we used the base including 
the states corresponding to much larger numbers of photons than for the states involved in 
\'^)cuf shall show that the analytical formulae (0) describing the evolution of the system 
restricted to the three states and the numerical results (with extended base) are in good 
agreement. This will ensure us that the applied method of the state truncation {non-linear 
quantum scissors) gives correct results. To prove it, we can analyse the fidelity defined by: 

F = \\{m\m)cut\\' , (8) 

where \4'(t))^^^ defined by (H)) is the truncated state and li'it)) is the state obtained by the 
numerical calculations for a large two-mode Hilbert space. Inspection of the time-evolution 
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FIG. 1: The time-evolution of the parameter 1 — F. The couphng strengths a = 7r/25, e = vr/25, 
and the nonhnearities Xa = Xb = 25. Time is scaled in 1/x units. 

of F enables us to estimate the quality of state truncation. The fidelity for the perfect 
truncation should give one. In particular, Fig.l shows the time evolution of the quantity 
1 — F. We can see that for the case discussed here, the deviations of F from the unity are 
only about 6 x 10""^, so the optical state truncation is fully justified. In consequence, our 
solutions indicate that for mode a the evolution is restricted to states corresponding to the 
vacuum, 1 or 2 photons, whereas for b the system evolves between the two states |0)^ and |2)j, 
only. Finally, we can get, to a high accuracy, the state (j3)), defined in the finite-dimensional 
Hilbert space, similarly as in P|. Therefore, the system described here can be referred to as 
two-mode nonlinear quantum scissors (TMNQS). Moreover, from the quantum information 
theory point of view, our model can be treated as a qutrit-qubit system ({0, 1, 2} ^ {0, 2}). 

III. ENTANGLEMENT AND THE BELL-LIKE STATES GENERATION 

In the previous section we have shown that the system's dynamics can be closed within 
a finite set of states, and hence, can be treated as TMNQS. This section is devoted to 
another very interesting and important feature, namely the system's ability to generate the 
maximally entangled sates (MES). 

Since the system evolution can be described in terms of the qutrit-qubit system, one can 
expect that for the model discussed here the Bell-like states can be generated. To investigate 
this phenomenon we plot the probabilities for the states involved in our model. Thus, Fig. 2 
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shows that these probabihties oscillate and some of the plots intersect for the values close 
to 0.5. This indicates that the Bell-hke states could be generated in the situation discussed 
here. In particular, we observe this crossing for the probabilities corresponding to the states 
|2)^ |0)j and |0)^ |2)^. Therefore, we shall concentrate on the states: 

l^i) = ^(|2)JO), + i|0)j2)J , 

\B,) = i=(|2)j0),-z|0)j2),) . (9) 

Moreover, we see that the probabilities corresponding to |2)j, and |0)^ |2)j exhibit the 
same behaviour (crossing for the values close to 0.5). This corresponds to generation of the 
states I Pi) and IP2), defined as follows: 

\Pi) = ^(Il)a|2), + ^|0)J2),) , 

1^2) = -^(|l)j2),-i|0)j2)J . (10) 

It is seen that these states can be expressed as 

|A,2) = ^(|l)±i|0))„®|2), , (11) 

where we get the product of two states, defined within various subspaces corresponding 
to modes a and b, respectively. In consequence, we cannot expect entanglement when the 
system reaches the states |Pi,2) and the product state is generated. 

To find whether the system can produce the states we are interested in, we plot (Fig. 3) the 
fidelities corresponding to the states I-B1.2) and |Pi,2)) ^-e. | {ip{t)\Bi^2) \ and | {ip{t)\Pi^2) \ 
respectively. We see, that for the properly chosen moments of time, these states can be 
generated to a high accuracy and some of the maxima in Fig. 3a are very close to unity. 
Moreover, we see that the fidehties corresponding to the Bell-hke states reach other maximal 
values, which are however, close to ~ 0.7. This indicates the existence of entanglement in the 
system, however, we get the states that are not 100% entangled. However if we enhance the 
coupling between the oscillators, the probability of creation of the state |2)^ decreases 
significantly and as a consequence, the fidelities for formation of the states |-Bi,2) are almost 
equal to unity (Fig. (3b) - solid line). Moreover, the oscillations between the states become 
faster and the system tends to be a qubit-qubit one. 

It is obvious that we could inspect other combinations of the pure states discussed here. 
For instance, the Bell-hke states involving |1)^ |2)j^ and |2)^ |0)j, could be considered. After 
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FIG. 2: The probabilities for the states |2)^ |0)^ (sohd hne), |0)^ |2)^ (dashed-dotted hne) and 
|2)j, (dotted hne). The couphng strengths a = 7r/25, e = 7r/25, and the nonhnearities Xa = 
Xb = 25. Time is scaled in 1/x units. 
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FIG. 3: The fidelities corresponding to the Bell-like states \Bi) (sohd hne), |i?2) (dashed line), and 
the product states (dotted line), IP2) (dashed-dotted line). The parameters at a) are identical 
to those of Fig.2. b) differ in e = 7r/5. 

inspection of the results plotted in Fig.2, we see that these Bell-like states could play some 
role in the system's evolution. We observe the crossings of the plots of the appropriate 
probabilities, however, their points of intersection correspond to the value of the probability 
equal to ~ 0.1. Therefore, for these moments of time the system is practically in the 
pure state |0)^ 12)^^. Nevertheless, one can see that for some other moments of time these 
probabilities become close to 1/2, although they do not intersect. Consequently, one can 
expect Bell-like states generation again. Fig. 4 shows the fidelities corresponding to the 
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FIG. 4: The fidelities corresponding to the Bell-like states: IB^) (dashed-dotted line), {B^} (dashed 
line), l^s) (solid line), \Bq) (dotted line). All parameters are identical to those of Fig.2. 

following Bell-like states: 
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We see that for the time ~ 12 the fidelity corresponding to the state {B^) becomes almost 
equal to unity. Hence, we get the Bell-like state IB^) slightly perturbed by the pure state 
|2)J0),. 

The ability for generation of maximally entangled states can be described in terms of 
the von Neumann entropy E of the reduced density matrix pa — Tri,pab or pi, = TraPab 
constructed on the basis of the full density matrix pah = \'ip) corresponding to the evolu- 
tion of the whole system. Equivalently, we can use the Shannon entropy H for the squared 
Schmidt coefficients Pi 

E = -TrpalogaPa = -TrpblogsPb 

^ H{pi} ^ -^Pilog^Pi . (13) 

i 

The entropy changes its value and is equal to zero for separable states, whereas for the 
maximally entangled states it should be equal to 1 ebit. As we see from Fig. 5, the entropy 
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reaches a limit of 1 ebit for several times for which the Bell-like (or close to them) states 
are formed. So, it proves that the NC with nonlinear interactions between field modes can 
be treated as a source of maximally entangled states. 




FIG. 5: The entropy of entanglement, in ebits, for the same parameters as in the previous figures. 



IV. EVOLUTION OF NONLOCALITY 

Both the entanglement of quantum states and the quantum nonlocality are essential 
and intriguing aspects of quantum mechanics. Quantum nonlocality manifests itself by a 
violation of the Bell- type inequalities Q]- It has been demonstrated that states which 



states. There exist also states (known as Werner ones 2l|) which are entangled nevertheless 



violate the Bell-type inequalities are the entangled pure states or a mixture of Bell 



they do not violate any of the Bell inequalities. It is also possible to obtain quantum states 
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for which we achieve the nonlocality without any entanglement 

Another problem is the way we can study the quantum nonlocality and its evolution. For 
this purpose we have usually used various specific parameters obtained numerically. For a 
two-qubit svstem we can examine the quantum nonlocality analysing the Bell-CHSH type 
inequality 

\Tr{pBcHSH)\<2 , (14) 
where Bell operator Bchsh has a form: 

BcHSH = a ■ a ® (b + b') ■ a + a' ■ or ® (b - b') • a . (15) 
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a, a', b, b' are the unit vectors in 3?^ and a is a vector of Pauli spin matrices cTi {i = 1, 2, 3). 
It is known that any state p can be expressed as: 



P=M/(8)/ + r- 0-(8)/ + /(g)S-0-+ ^ tnm<^n ® O^m | , (16) 
\ n,m=l / 

where r, s are also the unit vectors in 3?^, scalar product of r ■ a should be understood as 
Yli=i''^i^i ^'^d the coefficients tnm = Tr{pan ® crm) form a real matrix Tp. If we also define 
the real symmetric matrix Up = TjTp (Tj is the transposition of matrix Tp) we can use for 
further considerations the Miranowicz measure of nonlocality [2^: 

B{p) = vWr{MZGoy^T} . (17) 



M(p) as a measure of violation of inequality ((13)) has been proposed in j2J, |23] and takes 
the following form: M{p) = maxj^k{uj + u^}, where Ui {i = 1,2,3) are the eigenvalues of 
matrix Up. The Bell-CHSH inequality is then violated by the density matrix p if M{p) > 1. 
One can also study other parameters, for instance n{p) = max{0, M(p) — 1}, proposed in 
2^. The parameter B{p) defined in (jl7|) varies from zero to unity and is equal to unity for 
the case of maximal violation of the inequality (|T^ . and its value corresponds to the degree 
of Bell-CHSH inequality violation. As our model can be treated as a qutrit-qubit system we 
shall use I27I a local projection operation Ho,2 ® Ho,2 = (|0) (0| + |2) (2|) ® (|0) (0| + |2) (2|) 
to obtain the density operator for the qubit-qubit system. As a consequence, we get the 
reduced density matrix p^^^ of the form: 

^ ~ Tr (Ho,2 ® no,2pHo,2 ® Ho,2) ' ^ ^ 

Here, we are able now to analyse the degree of quantum nonlocality using (jl7|) for the two 
qubit system obtained from two qutrit system by the above projection. As an example, 
we show the time dependence of the B{p) applied for p^^^ and trace its value for the times 
where Bell-like states (eq.®) are formed. Comparing this figure with fig. 3a we can easily see 
that the maximal degree of Bell-CSHS inequality violation appears when one of the Bell-like 
states (eq.Q) is generated with perfect fidelity. The maximal violation also appears when 
one of the Bell-like states is formated with less fidelity (maxima labelled as 2° and 3° at the 
fig. (6a)). Maximum 2° for example appears when the fidelity of formation of |i?2) state is 
zero but at the same time there appears a state \Pi) with large fidelity - however, \Pi) has 
also the |0)^ |2)^ component. Thus, we would say that the infiuence of the states specified in 
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FIG. 6: The Bell inequality violation degree corresponding to the formation of Bell-like states 
\Bi^2) for the coupling strengths e = 7r/25 (a) and e = 7r/5 (b). Time is scaled in 1/x units. The 
labels {1*^, 2*^, 3*^} enumerate the maxima - see the discussion in the text. 

eq.(|l(jp give rise to that nonlocality effects. Larger values of the coupling strength result in 
faster oscillations of the parameter B{p) considered, similarly as for the fidelities shown in 
Fig. 3b. In consequence, the Bell-CSHS inequality is maximally violated when the Bell-like 
states \Bi^2) form (Fig.(6b)). 

V. THE DAMPING CASE 

As mentioned earlier, the damping processes are able to modify physical properties of 
systems containing nonlinear oscillators jlS^. In particular, the ability of finite dimensional 
states formation could be considerably reduced. Therefore, we shall concentrate in this 
section on the influence of damping processes on the Bell-like states generation, and we will 
inspect the effect of the damping process (caused mostly by the losses at cavity mirrors) on 
the formation of maximally entangled states. Therefore, we will have to include into the 
description of the system dynamics the appropriate quantities corresponding to the energy 
leakage from the cavity. It is obvious that this problem concerns both a and h modes. 
Consequently, we have to define two collapse operators corresponding to the energy losses 
in these two modes. They are: 



C 



1 — 





(19) 
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(21) 



We will now turn to the standard methods of the quantum theory of damping, to find the 
evolution of the damped system by numerically solving the master equation for the density 
matrix. The generic form of the master equation is 

(20, 

where the Liouvillian £ is a superoperator and in our case it takes the form: 
t = -{Hp-pH)+Y, CkpCl - - (clCkP - pClCu 

k=l L 

Once we have found the Liouvillian, we are able to solve the master equation and get the 
time evolution of the density matrix elements, which corresponds to the states involved 
in the dynamics of the coupler described. The evolution is found as a series of complex 
exponentials exp(sjt), where each Sj is an eigenvalue of C. 

Although it is possible to examine the effect of the damping processes on the generation of 
all the entangled states discussed here, we shall concentrate in this section on the generation 
of and |i?2) states. We assume that the processes discussed here perturb the mechanisms 
of other Bell-like states generation in the same way. Moreover, we shall compare our results 
with those concerning NC with a linear internal coupling Q| . Therefore, we have assumed 
that the values of the parameters describing our system are comparable with those discussed 
in - the values of the parameters discussed there are assumed to be more realistic from 
the experimental point of view. Since this part of the paper is devoted to the damped case 
and we use the description of the system based on the density matrix formalism, we use the 
following definition of the fidelity F: 

F = TrVp'/'Pcutp'/' , (22) 

where p and pcut are the density matrices corresponding to the states compared. For the 
case discussed in the previous section this fidelity becomes identical with that defined in 
eq.®. 

Thus, Fig. 7 shows the fidelities corresponding to the states \Bi) and |i?2) for different 
values of the damping parameters. We see that for a very weak damping case we get Bell- 
like state |i?2) almost perfectly (the same situation can be observed for the model discussed 
in ^^). However, if we increase the values of the damping parameters the effectiveness 
of the process becomes considerably reduced. Nevertheless, we see that for k = we 
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FIG. 7: Fidelities corresponding to the Bell-like states \Bi) (solid lines) and |i?2) (dashed lines) for 
various values of the damping constants Ka = Kb = k: k = x/500 (a), k = x/75 (b), and k = 
(c). The remaining parameters are: Xa = Xb = X = lO^rad/sec, a = x/20, e = a/2. 

get the desired state quite effectively. One should note that for this case, not only the 
maximal value of the fidelity corresponding to the state I-B2) is important, but also the rate 
of simultaneous vanishing of the fidelity for the state \Bi) should be considered, due to the 
fact that the sum of these states gives the pure state |2)^ \0)^. For the case of k = x/75 
both conditions are sufficiently fulfilled . The strengths of the damping discussed here are 
effectively greater than in and consequently, our model seems to be more suitable for the 
Bell-like states generation, as it is less sensitive to the losses in the system. We can explain 
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this fact as a result of the resonant nature of the nonhnear internal coupling. The coupling 
with the external reservoir has a linear character and involves the one-photon and vacuum 
states corresponding to two modes a and b, and then the states discussed here ({|0)^ |2)^, 
|2)^ |0)^, |1)^|2)^}). Consequently, due to the resonant character of the transitions inside 
the set of these three states, the transitions become dominant and couplings to the states 
interacting with the external bath can be neglected for the cases considered here. This result 
is a consequence of the nature of our system that is constructed on the nonlinear scissors 
model base, where the resonances play a crucial role Q]. Once the damping process causes 
the decay of the fidelity for the Bell-states, one can think about the possibility of applying 
apurification protocol to reduce the influence of the damping. It has been already shown 
22 . [2^ that there is no local purification scheme which produces a maximally entangled 



state from a mixed state entangling two finite-dimensional systems. Therefore, one should 
concentrate on the global purification techniques. For these cases the model discussed here 
would be only a part of a greater system, in which the information would be shared between 
two parties and measured by them. Moreover, one should keep in mind that practically 
all systems that are discussed in that context in literature are two-qubit systems, while 
the ours is a qutrit-qutrit one. However, Alber et.al 29| have discussed higher dimensional 
systems, and this fact gives us some hope for finding the appropriate purification protocol 
that would help improve the maximally entangled state generation mechanism. They have 
proposed an application of generalized XOR-gate that can be constructed on the basis of 
nonlinear Kerr interaction governed by the Hamiltonian H = xnin2 (that differs from those 
discussed in our paper). In the context of the purification protocol presented there, the 
model discussed here can be treated as promising proposal for further investigation of the 
methods of maximally entangled states generation. The system discussed here could be a 
part of a more complicated model that allows generation of maximally entangled states with 
the help of the above mentioned purification method. 



VI. CONCLUSIONS 



We have discussed the model of a coherently driven NC with a nonlinear internal coupling 
showing that for sufficiently small values of the coupling constants we can apply the nonlinear 
scissors method. Consequently, the system's evolution becomes closed inside a set of three 
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states. Moreover, which seems to be the most interesting, we have shown that our system 
can be treated as a maximally entangled states generator, in particular, it can generate the 
Bell-hke states. We have also discussed the fragility of the system to the damping processes. 
The model proposed here, due to the nonlinear character of the internal coupling, seems 
to be more suitable for practical realization than that involving only linear interactions. 
As we discuss the context of practical feasibility of our model we should keep in mind 
the assumptions of very weak interactions. Consequently, the model could be practically 
realizable if we assume sufficiently large values of the nonlinearities discussed here. This 
fact indicates that the most suitable for practical realization would be systems based on the 
electromagnetically induced transparency phenomenon or atomic dark resonances. These 



have been proposed by Schmidt and Imamoglu 
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32| | and experimentally 



processes 

observed [33[. For models involving these effects one could expect the possibility of producing 
the so-called giant Kerr nonlinearities. Imamoglu has estimated that for such systems 
we can achieve x ~ 10*^ rad/sec 0, 1^^. So, we see that the system discussed in this 
paper can be treated at least as a potential source of the Bell-like states. Moreover, it 
could be implemented as a part of a more complicated model interesting from the quantum 
information theory point of view, keeping in mind all obstacles and limitations concerning 
its practical realization. 
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